Casimir spring and compass: Stable levitation and alignment of compact objects 
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We investigate a stable Casimir force configuration consisting of an object contained inside a 
spherical or spheroidal cavity filled with a dielectric medium. The spring constant for displacements 
from the center of the cavity and the dependence of the energy on the relative orientations of the 
inner object and the cavity walls are computed. We find that the stability of the force equilibrium 
can be predicted based on the sign of the force, but the torque cannot be. 



The Casimir force between atoms or macroscopic ob- 
jects arises from quantum fluctuations of the electrody- 
namic field [H [2] . In all known examples it is found that 
the force is attractive, as long as the space between the 
objects is empty, and the magnetic susceptibility of the 
objects is negligible compared to their electric suscep- 
tibilities. But when space is filled with a medium with 
electric permittivity e m intermediate between that of two 
objects, ei < cm < £2, the force between the two becomes 
repulsive j3] . This effect has recently been verified exper- 
imentally |3]. But while repulsive forces are nothing new, 
they become interesting for applications when they pro- 
duce stable equilibria, which, for example, the Coulomb 
force cannot. 

For an infinite cylinder enclosed in another, the 
Casimir force has recently been shown to have a stable 
equilibrium in the two directions perpendicular to the 
cylinder axes, when the material properties are chosen 
so that the force between two slabs of the same materials 
would be repulsive [5]. If the inner and outer cylinders 
have square cross sections, the direction of the torque 
exerted by one cylinder on the other is found to agree 
qualitatively with the predictions of the pairwisc summa- 
tion or proximity force approximations (PSA or PFA). 
Orientation dependence of the Casimir energy has also 
been studied recently for a small spheroid facing another 
spheroid or an infinite plate [6], or for a small spheroid lo- 
cated inside a spherical metal cavity [7] . (The smallness 
of the spheroid refers to keeping only the first term in 
the series expansion of the Casimir energy in the largest 
length scale of the spheroid.) The preferred orientation 
of the small spheroid flips when it is moved between the 
inside and the outside of the spherical shell; there is no 
torque if the spheroid faces an infinite metal plate, which 
is the limit of an infinitely large shell. 

In this paper we investigate the first configuration, de- 
picted in Fig. [T] in which a compact object levitates sta- 
bly due to the Casimir force alone. We characterize the 
Casimir energy 

£ = £ Q + \k^ + ^k^ + ^k^ + --- (l) 

by the spring constant k and the coefficients k n in a se- 
ries expansion in a/R, where a is the magnitude of the 
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FIG. 1: (Color online) Summary of the configurations we con- 
sider and of the results. To be concrete, we have assumed 
that the small spheroid's zero frequency permittivity satisfies 
ei.o > £m,o and that it is larger in the body-fixed z direction, 
so af z > j_. Furthermore, the magnetic permeabilities are 
all set to one. a) Direction of the force on such a spheroid in 
a spherical cavity if f\ > 0, which holds when £m,o > eo,o, 
and the direction of the torque when f% > 0, which is satis- 
fied when either £a/,o > eo,o or eu,o <S eo,o- b) A finite size 
sphere experiences a restoring force for the various combina- 
tions of materials listed in Table [I] c) Direction of the torque 
in the center of a slightly spheroidal cavity if g E > 0, which 
requires either e M ,o < eo,o or e M ,o > eo,o- 

displacement from the center of the cavity and R the ra- 
dius of the (undeformed) spherical cavity, (k has units of 
energy here.) Unlike the case of infinite cylinders, where 
one object can be displaced along the cylinder axes with- 
out changing the energy, our case exhibits, for appropri- 
ately chosen materials, true stability in all directions and 
applies to realistic situations. For example, we compute 
the force on a metal sphere in a spherical drop of liquid 
surrounded by air. By measuring the mean square de- 
viation from the center, (a 2 ) — ikgTR 2 /k, the spring 
constant k can be verified experimentally. We can es- 
timate that the size of the droplet R has to be smaller 
than w 3/im-^3 , where r is the typical dimension of the 
inner object, for the thermal motion to be confined near 
the center of the cavity. (This length scale is obtained by 
balancing SksT for room temperature with a rough esti- 
mate of the spring constant, k ~ ~ ■) To keep the two 
objects nearly concentric against the gravitational force, 
R has to be smaller than ss 1/im for the typical metal or 
liquid densities considered here. A variety of applications 
may benefit from our analysis; in cancer therapy new 
treatments utilize nanocarriers that trap drug particles 
inside R ~ .1/im polymer or lipid membrane shells [H], 
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and molecular cages are proposed as containers for the 
storage of explosive chemicals [9] . Our results may guide 
the search for better materials and sizes of the enclosing 
cell. 

Whether the center of the cavity is a point of stable 
or unstable equilibrium turns out to depend on whether 
the Casimir force is repulsive or attractive for two sep- 
arated objects under those conditions. The direction of 
the torque, on the other hand, depends on the dielectric 
properties of the medium and the cavity walls in an un- 
intuitive way, which cannot be predicted by the PSA or 
PFA, see Fig. [l] In particular, this behavior is not due 
to dispersion effects, which explain similar phenomena 
reported in Ref. [5]. We calculate the torque on a small 
spheroid that is displaced from the center of a spherical 
cavity (Fig. Ilia)) or concentric with a slightly spheroidal 
cavity (Fig. IT] c)). In the former case, the orientation 
dependence manifests itself in k. In the latter case, £q 
depends on the relative orientation of the spheroid in- 
side the deformed cavity. (If the cavity is spherical, £q 
is a finite constant that can be ignored.) By choosing 
appropriate materials, the inside object, e.g., a nanorod, 
can be made to align in different ways with the cavity 
shape, a situation which is reminiscent of a compass nee- 
dle aligning with the magnetic field of the earth. 

The starting point of the analysis is the scattering the- 
ory approach. The method is explained and derived in 
detail in Ref. |10j , where a partial overview over its pre- 
cursors, e.g. [TTJ rjJl H3] , is provided. The Casimir energy 

he f°° ■ 
£= — d/dndet(T-^W°JFf e V 10 ) (2) 
2ti" Jo 

is expressed in terms of the inner object's exterior T- 
matrix, J-j e , and the outer object's interior T-matrix, 
J-q. The exterior T-matrix describes the scattering of 
regular wave functions to outgoing waves when the source 
lies at infinity. The interior T-matrix expresses the oppo- 
site, the amplitudes of the regular wave functions, which 
result from scattering outgoing waves from a source in- 
side the object. The translation matrices W w and V l ° 
convert regular wave functions between the origins of the 
outer and the inner objects; they are related by complex 
transpose up to multiplication by (—1) of some matrix 
elements. 

With uniform, isotropic, and frequency-dependent per- 
mittivity e x (icK) and permeability fi x (icn) functions (x — 
I: inner object; x — O: outer object; x = M: medium) 
the T-matrix of the sphere is diagonal. It is given by 

for E (electric) polarization and by the same expres- 
sion with ^p- replaced by — for M (magnetic) po- 



larization (not to be confused with subscript M indi- 
cating the medium's response functions). In Eq. ^ 
the frequency dependence of the response functions has 
been suppressed. The indices of refraction n x (icK) — 
■>J e x (icn)ii x (icK) of the sphere and the medium appear 
in the ratio zi(ick) — ni(icn) / 'tim^ck) and the argument 
£ = riM(icK)Kr. The first equality in Eq. ^ defines an 
abbreviation for the T-matrix, in which the superfluous 
polarization and angular momentum (/, m) indices are 
suppressed. The interior T-matrix of the spherical cavity 
is obtained from the exterior T-matrix of the sphere by 
inserting the outside object's radius and response func- 
tion in place of those of the inside object and exchanging 
the modified spherical Bessel functions ii and ki every- 
where. 

However, the scattering approach is not limited to 
simple geometries. An array of techniques is avail- 
able for calculating the scattering amplitudes of other 
shapes. We employ the perturbation approach to find 
the T-matrix of a deformed spherical cavity [21 [TS] of 
radius R + 6(1 — 3/2sin 2 #). The deformation, indi- 
cated in Fig. [l]c), is chosen so that the volume is un- 
changed to first order in d. We find the O(S) correction, 
to the T-matrix in a perturbation series expansion, 
T = J-^ +J-^ + ■ ■ • , by matching the regular and outgo- 
ing fields according to the Maxwell boundary conditions 
along the deformed object's surface [IB]. On the other 
hand, for a small object (compared to the wavelength 
of the radiation), we can approximate the T-matrix to 
lowest order in k using the static polarizability tensor, 
^iAmP.im'P = 2/3(n M:0 K) 3 a^ m ,+O(K 5 ), where the sub- 
script indicates the static (icn — 0) limit and P is 
the polarization label. The T-matrix elements involving 
higher angular momenta I > 1 are higher order in n. For 
a small ellipsoid, in particular, the electric polarizability 
tensor ot E is diagonal in a coordinate system aligned with 
the ellipsoid's body axes, 

E V £ I,0 ~ e Af,0 ,.s 

a u = -, r — , (4) 

47T e M<0 + (e/,0 - £A/,o)«-i 

where i € {x, y, z} |17j . The larger the semi-axis in di- 
rection i, the smaller the depolarization factor rii (not to 
be confused with the index of refraction). The entries 
of the magnetic polarizability tensor a M are obtained 
by exchanging /u Xj o for e x $ in Eq. Q. In the small size 
limit the polarizability tensor of a perfect metal ellipsoid 
is obtained by taking both e^o to infinity and /z/.o to 
zero. 

For simplicity we specialize to a spheroid, which has 
two equal semiaxes. We choose the semiaxes along x 
and y to be equal, therefore, a xx — (Xy y = a^ ± . The 
direction of displacement of the spheroid from the center 
of the cavity is the lab's z axis. 9 denotes the angle 
between the spheroid's and the lab's z axes. For such 
a small spheroid inside a spherical cavity of radius R, 
the spring constant is obtained by expanding the log- 
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determinant in Eq. ([2| to first order, 
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express the rotation invariant and the orientation depen- 
dent parts of the energy, respectively, f^ 1 and f^ 1 are 
obtained from ff and jvf by exchanging E and M ev- 
erywhere. This result is valid for asymptotically large 
R; it involves only the zero frequency (icn = 0) response 
functions and the 1 — 1,2 scattering amplitudes of the 
cavity walls in Eqs. ([5| and ([6]). 

The behavior of the functions ff, depicted in Fig.[2]for 
Mm,o == Mo,o, is as expected: ff is monotonic, positive 
when 6m,o > eo,o, negative when ejy-,0 < e o,o, and /j^ 
always has the opposite sign of ff. When /-p is positive, 
a small object with > ^m,o is levitated stably, when 
/j 6 is negative, cm,o > e /,o has to hold. Thus, stability 
occurs under the same conditions as repulsion for two ob- 
jects outside of one another. The opposite sign of f^ 1 is 
expected from equivalent expressions for the two-infinitc- 
slab geometry [18]. When the dielectric contrast between 
the medium and the outer sphere is taken to small or 
large limits, stability or instability is maximized. To ver- 



filled with various liquids. [23J The coefficients k, k^, and 
fcg in the series expansion in Eq. ([IJ are listed in Table [I] 
For comparison, the asymptotic result kn^oc is also in- 
cluded. If both inside and outside object are spherically 
symmetric, the series in Eq. ([!]) does not contain terms 
with n odd. 
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nations of materials for the case depicted in Fig. [T]b). The 
dimensionless numbers in the table have to be multiplied by 
R is given in microns [/im], fc.R->oo depends on R only 
through the ratios ^ and ^ , so its numerical prefactor is the 
same for all 7?. The highest cutoff used was i max = 30. 
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FIG. 2: (Color online) ff and ff^ describe the part of the 
spring constant ku^^o, which is invariant under a rotation 
of the inside object. The vertical lines indicate the values 
pertaining to the configurations presented in Table [I] ethanol- 
vacuum (0.16), bromobenzene- vacuum (0.30), and gold cavity 
walls (1). In this plot, ^ia/,o = ^o,o- 

ify whether stability is observable for realistic materials 
and object sizes, we evaluate Eq. ^ numerically for a 
sphere of radius r inside a spherical cavity of radius R 



The three materials were chosen so that the sequence 
of permittivities e/, cm, £o cither increases or decreases 
for the imaginary frequencies that dominate in Eq. Q. 
Contrary to the prediction of the PSA and PFA the force 
is not symmetric with respect to exchange of the inner 
and outer permittivities. In the same medium, a high di- 
electric sphere is held more stably in the center of a cav- 
ity with low dielectric walls than a low dielectric sphere 
inside a cavity with high dielectric walls. 

The asymptotic result kji^^ yields a good approxi- 
mation of k for = \- But from Eq. (p| one would 
expect fc_R^oo to grow linearly with the volume of the 
inner sphere, since polarizability is proportional to the 
volume, see Eq. In fact, k for ~ = | is about 1000 
times larger than for = j, instead of just 27 times. 
This means that for a gold sphere in a liquid drop with 

— | and R = 1/im at room temperature, indeed, the 
Casimir spring holds the particle near the center effec- 
tively, y/fa?)/R < 0.1. 

Compared to the stability conditions studied thus far, 
the orientation dependence of the energy is more varied. 

and ff 1 , plotted in Fig. [3] have the same sign for 
most ratios of medium to outside permittivities, unlike 
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fi and fi 1 , which always have opposite signs. In these 
ranges of values, the contributions to the torque from 
electric and magnetic polarizability are opposite for a 
small perfect metal spheroid, for which e/ > ea/.o an d 
M/,o < Mm,o- 



Unlike /f and f™, also, fi and / 2 M change 



sign again at 



e M,0 



80 and at ^ 

e M,0 



2000, respectively. 



So, while the direction of the total force and the stability 
of the equilibrium can be predicted based on the relative 
magnitudes of the permittivities, the torque cannot be. 
The second sign change of and does not agree 
with the PSA or PFA. Furthermore, it raises the question 
whether calculations of the Casimir torque for infinite 
conductivity metals are 'universal' in the sense that they 
produce the correct qualitative results for real materials. 




FIG. 3: (Color online) f% and /«j describe the part of the 
spring constant kn->oo, which changes with the orientation of 
the inside spheroid. In this plot, (mm,o = HO,o- 

The orientation dependent part of the energy for the 
configurations discussed thus far vanishes, of course, 
when the small spheroid is in the center of the spheri- 
cal cavity. If the cavity is slightly deformed, however, 
the energy, £q, depends on the relative orientation of the 
spheroid and the cavity. We deform the spherical cavity 
as described earlier and obtain to first order in 5/R, 
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£ — > M] , (7) 



where the orientation-independent part of the energy has 
been dropped. For /j,o,o = Ma/,o, g E and g M are given by 
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where the arguments of the modified spherical Bessel 
functions and of zo(0) are suppressed, ii and ki are func- 
tions of £, and ki stands for ki(zo(0)£), where zo{icn) = 
no(icK)/riM(icK) is the ratio of the permittivities of the 
cavity walls and the medium. The material dependent 
functions g E and g M are plotted in Fig. |3J 




FIG. 4: (Color online) g E and g M describe the dependence of 
the energy So on the relative orientation of the inside spheroid 
and the deformed cavity walls. 



Again, g E has the same sign for 

Fig.H and 
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(left in 
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oo (right). In addition to the root at 
0.46. 



— 1, q E also vanishes at 

The rich orientation dependence of the energy is ex- 
pected to collapse as the size of the inside object grows 
to fill the cavity and the PFA becomes applicable. Based 
on the stability analysis for finite inside spheres, though, 
we expect the asymptotic results, f£ and g p , to predict 
the orientation dependence for reasonably small inside 
spheroids. 

In reality, various corrections to the idealized shapes 
have to be taken into account. A drop of liquid, which is 
placed on a surface, is influenced by gravity and interac- 
tions with the substrate. An asymmetric inner object, or 
one that is displaced from the center, deforms the shape 
of the droplet additionally by causing uneven Casimir 
stresses. 
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